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ABSTRACT 

A  limited  optimization  study  is  undertaken  of  imperfect,  eccentrically 
stiffened  2024-T3  aluminum  alloy  cylinders  under  axial  compression  based 
on  (l)  a  linear  theory  with  applied  buckling  load  reduction  factors  to 
account  for  initial  imperfections  and  inelasticity  and  (2)  a  kinematically 
and  constitutively  nonlinear  maximum  strength  analysis.  The  latter  is 
shown  to  be  less  conservative  for  the  most  efficient  design:  combined 
stringer-  and  ring-stiffening.  Credibility  of  the  maximum  strength  ap¬ 
proach  is  demonstrated  through  its  use  in  correlating  experimental  re¬ 
sults  obtained  for  10-  and  lk-in.  diameter  steel  cylinders  having  differ¬ 
ent  proportional-limit  stresses,  respectively.  Finally,  material  non¬ 
linearity  is  shown  to  offer  a  plausible  qualitative  explanation  for 
a  significant  amount  of  the  scatter  in  axial  compression  tests  for 
cylinders  with  radius-to-thickness  ratios  of  practical  interest. 
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FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  effort  being  undertaken 
at  Stanford  University  for  the  Eustis  Directorate,  U.  S.  Army  Air  Mobility 
Research  and  Development  Laboratory  under  Contract  DAAJ02-70-C-0075  (Task 
1F061102AJ3F02).  The  present  program  is  one  in  a  series  aimed  at  es¬ 
tablishing  accurate  theoretical  prediction  capability  for  the  static  and 
dynamic  behavior  of  aircraft  structural  components  using  both  conventional 
and  unconventional  materials.  Predecessor  contracts  supported  investiga¬ 
tions  which  led,  in  part,  to  the  results  presented  in  Reference  10,  11, 

23,  36,  40,  42,  and  43. 

For  his  part  in  assisting  with  digital  computer  calculations,  the  efforts 
of  Mr.  Joseph  Mullen,  Jr.,  are  gratefully  acknowledged. 
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INTRODUCTION 

The  post-World  War  II  developments  in  aircraft  design  leading  to  greatly 
increased  performance  and  capability  can  be  attributed  more  to  break¬ 
throughs  in  the  areas  of  propulsion  and  aerodynamics  than  to  advances 
in  structures.  With  weight-consciousness  becoming  ever  more  prevalent, 
it  is  interesting  that  the  primary  source  of  vehicle  empty  weight,  the 
load-carrying  structure,  continues  to  evidence  prejet-age  design  con¬ 
siderations.  Much  optimism  has  been  placed  in  the  potential  of  com¬ 
posite  materials  for  improving  structural  efficiency.  However,  the 
promise  of  composites  notwithstanding,  conventional  sheet-stringer 
construction  continues  to  form  the  basic  structure  of  fuselages,  wings, 
and  empennages  in  general;  thus,  it  offers  a  logical  and  attractive 
area  for  fundamental  weight  reduction,  with  possible  cascading  effects 
relative  to  the  overall  system  either  in  decreased  weight  or  weight 
tradeoffs. 

In  conventional  aircraft  construction,  the  thin  curved  plate  or  shell 
is  the  key  element.  Yet  the  design  of  such  elements,  especially  for 
stiffness  considerations,  continues  to  be  founded  upon  semiempirical 
and,  apparently,  overconservative  criteria.  The  buckling  of  thin 
shells  in  compression  is  a  case  in  point.  The  unfortunate  lack  of 
experimental  data  to  guide  the  development  of  satisfactory  design  pro¬ 
cedures  for  practical  shells  is  emphasized  by  Peterson  in  Reference  1. 
Additionally,  the  frustration  underlying  the  myriad  attempts  to  resolve 
the  thin  shell  buckling  problem  satisfactorily  from  a  theoretical  stand¬ 
point  is  evidenced  in  the  remarks  of  Almroth  et  al .  in  Reference  2  as 
justification  for  the  need  to  establish  approximate  methods  of  stability 
analysis  leading  to  design  criteria  which  are  less  conservative  than 
those  in  common  usage  but  which  can  still  be  employed  with  a  high  degree 
of  confidence. 

Much  has  been  made  of  the  deleterious  effects  of  initial  imperfections, 
boundary  conditions,  and  prebuckling  deformations  in  the  buckling  of 
shells  under  compressive  loading  in  attempts  to  explain  the  significant 
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discrepancy  between  the  60-year-old  classical  theory  and  experiments 

^  4 

of  the  past  50  years  (see,  for  example,  the  survey  papers  of  Hoff,'' 
Stein^,  and  Hutchinson  and  Koiter^).  The  consensus  leans  heavily 
toward  initial  imperfections  as  the  significant  degrading  factor 
compared  to  the  influences  of  boundary  conditions  and  prebuckling 
deformations.  The  importance  for  design  purposes  attributed  to  ini¬ 
tial  imperfections  in  the  buckling  of  thin-walled  shells  in  compression 
is  reflected  in  recommended  criteria  appearing  in  handbooks  either  de¬ 
veloped  or  sponsored  by  the  National  Aeronautics  and  Space  Administra- 
7  8  9 

tion.  ’  Also,  one  of  the  most  recent  of  design  criteria  advanced 
for  axially  loaded  cylindrical  shells  makes  generous  use  of  an  empir¬ 
ical  factor  to  account  for  the  presence  of  initial  imperfections  in 
2 

practical  shells. 

Although  from  time  to  time  the  question  of  nonlinear  elastic  and  off¬ 
set  plastic  (yielding)  buckling  has  arisen,  the  tendency  has  been  to 
treat  the  buckling  of  thin-walled  highly  curved  plates  and  shells  on 
the  basis  of  mathematical  models  using  linear  stress-strain  laws. 

Plastic  buckling  of  shells  is  more  appropriate  to  thick  shells  and 
of  little  interest  to  aircraft  shell  designers  and  analysts.  Very 
thin  shells,  wherein  linear  elastic  material  behavior  at  buckling  is 
more  prone  to  occur,  is  again  of  diminished  interest  in  design  because 
of  the  large  penalties  to  be  paid  for  imperfection-sensitivity  "knock¬ 
down"  factors  and  loss  of  aerodynamic  smoothness  due  to  the  catastrophic 
nature  of  thin  shell  buckling.  In  the  range  of  "effectively  thick* 
shells  (closely  stiffened  and  sandwich  construction),  however,  the  im¬ 
perfection  sensitivity  is  significantly  less.  But,  it  is  in  this  range 
that  practical  configurations  exist  (closely  spaced,  moderately  heavy 
stiffening  and  sandwich  concepts),  and  their  buckling  loads  cannot  be 
determined  independently  of  nonlinear  material  behavior  if  structural 
efficiency  is  to  be  achieved  subject  to  the  usual  condition  that  no 
yielding  take  place  at  limit  load.  Mayers  and  Wesenberg^'  ^  have 
called  attention  to  the  importance  of  nonlinear  material  behavior  in 
the  failure  mechanism  of  imperfect,  axially  compressed,  conventional, 
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stiffened,  and  sandwich  shells,  especially  in  the  range  of  practical 
interest  for  both  material  properties  and  geometric  parameters. 

It  is  the  purpose  of  the  present  study  to  demonstrate  how  the  maximum 
strength  analysis  can  be  used  to  optimize  in  a  broad  sense  an  initially 
imperfect,  eccentrically  stiffened,  aluminum  cylinder.  The  resuJ's  are 
compared  with  the  more  common  procedure  in  Reference  9*  Additionally, 
credibility  of  the  maximum  strength  concept,  theory,  and  analysis  is 
demonstrated  by  application  to  several  of  the  steel  test-cylinder  con¬ 
figurations  reported  in  Reference  12;  these  cases  serve  as  an  explanation 
of  the  apparently  anomalous  cylinder  behavior  appearing  in  the  presenta¬ 
tion  of  results  in  Reference  12  under  the  assumption  that  all  buckling 
failures  occurred  with  linear  elastic  material  behavior.  Finally,  the 
overall  implication  of  the  present  investigation  and  its  forerunners ^  ^ 

possibly  sheds  new  light  on  interpretation  of  the  results  of  the  per- 

13  14 

ceptive  pioneering  work  performed  by  Donnell  and  Lundquist  both 
theoretically  and  experimentally  on  the  compressed  shell  buckling 
problem  almost  40  years  ago  and,  apparently,  too  soon  overlooked  by 
myriad  investigators  concerned  with  mathematical  solutions  of  the 
classical  small-  and  large-deflection  theories  for  linear  elastic 
shells,  including  such  effects  as  initial  imperfections,  combinations 
of  boundary  conditions,  and  prebuckling  deformations. 
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DISCUSSION  AND  RESULTS 


TYPICAL  DESIGN  PROCEDURE 


The  buckling  phenomenon  attributed  to  shells  under  compressive  loading 
has  been  commonly  treated  in  analysis  and  design  as  classical  bifurca¬ 
tion  buckling  with  empirically  and  theoretically  established  reduction 
factors  applied  to  the  classical  buckling  loads  to  account  for  the  sig¬ 
nificant  discrepancies  that  invariably  exist  between  theory  and  experi¬ 
ence.  The  classical  axial  compression  buckling  loads  for  circular 
orthotropic  shells  with  closely  spaced,  eccentric  stiffening  can  be 
obtained,  for  example,  from  the  initial  developments  of  Baruch  and 
Singer*"^  and  the  broad  generalization  accomplished  by  Block  et  al.*® 

For  symmetric  stiffening,  the  Baruch-Singer  buckling  criterion  reduces 

17 

to  that  of  Stein  and  Mayers  for  orthotropic,  circularly  curved  sand¬ 
wich  plates  and  shells  when  transverse  shear  effects  are  neglected. 

The  latter,  in  turn,  again  with  transverse  shear  effects  neglected, 
reduces  to  the  classical  buckling  load  for  homogeneous,  isotropic, 
long  thin  shells  developed  independently  by  Lorenz,*®  Timoshenko, *^ 
and  Southwell  some  60  years  ago,  For  such  cylinders,  prediction  and 

experience  have  finally  been  closely  matched  by  very  careful  fabrica- 
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tion  and  testing  techniques.  ’  ’  ’  Although  the  specimens  in¬ 

volved  must  be  considered  as  much  too  small  for  practical  applications, 
they,  nevertheless,  served  to  expose  initial  imperfections  as  a  main 
contributor  to  the  discrepancy  between  classical  theory  and  experiment 
as  well  as  to  explain  the  cause  of  some  of  the  significant  scatter  ex¬ 
isting  in  a  very  large  body  of  test  data. 

Buckling  Load  "Knockdown"  Factors 

One  of  the  most  scrupulous  attempts  to  establish  a  consistent  "knock¬ 
down"  factor  curve  as  a  function  of  radius-to-thickness  ratio  for  use 
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in  design  is  evidenced  by  the  work  of  Weingarten  et  al .  The  curve 
for  axially  compressed  cylinders  presented  in  Reference  25  represents 
the  lower  bound  to  selected  buckling- load  data  on  unstiffened  shells 
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garnered  from  a  host  of  sources.  Available  data  were  eliminated  from 
consideration  when,  for  example,  the  test  results  were  believed  to  re¬ 
flect  the  influence  of  boundary  conditions,  poor  fabrication  and  testing 
techniques,  and  average  stress  levels  in  excess  of  70$  of  material  yield 
stress.  The  selection  process  notwithstanding,  the  group  of  data  leading 
to  the  lower  bound  design  curve  still  presents  a  considerable  scatter 
band  at  any  radius-to- thickness  ratio,  not  all  of  which  in  the  opinion 
of  Weingarten  et  al.  can  be  isolated  to  the  effects  of  initial  imper¬ 
fections.  The  design  curve  of  Reference  25  is  essentially  the  same  as 
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those  appearing  in  the  NASA  handbooks  ’  ’  for  use  with  unstiffened 
cylinders.  Careful  attempts  have  also  been  made  to  explain  the  sig¬ 
nificant  scatter  band  obtained  with  stiffened  cylinders  (see,  for  ex¬ 
ample,  the  study  by  Peterson  in  Reference  l).  After  such  considerations 
as  prebuckling  deformations,  discreteness  (versus  smearing)  of  stiff¬ 
ening  elements,  and  fabrication  tolerances,  the  prasence  of  imperfec¬ 
tions,  residual  stresses,  and  local  skin  buckling  are  the  leading  reasons 
advanced  for  the  lack  of  correlation  between  classical  theory  predictions 
and  experiments  relative  to  compressed  shells.  Unfortunately,  the  number 
of  stiffened-cylinder  tests  reported  in  the  literature  is  miniscule  com¬ 
pared  with  that  for  unstiffened  cylinders;  further,  as  pointed  out  by 
Peterson, 1  there  are  no  available  data  on  conventional  ring-  and  stringer- 
stiffened  cylinders  in  which  local  skin  buckling  does  not  precede  general 
buckling. 


Alternate  Procedure 
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In  the  case  of  either  stiffened  or  unstiffened  shells  under  compression, 

there  has  developed  more  recently  an  alternative  approach  to  that  of 

applying  suitable  reduction  factors  to  the  classical  buckling  load. 

This  approach  stems  from  the  Koiter  theory  for  the  elastic  buckling 

and  initial  postbuckling  of  imperfect  thin  shells.  In  Koiter-type 
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analyses,  ’  the  calculated  buckling  load  inherently  accounts  for 
the  presence  of  an  imperfection  distribution  over  the  shell  surface 
corresponding  to  the  axisymmetric  classical  buckling  mode.  The  imper¬ 
fection  distribution  possesses  an  arbitrary  amplitude  which  can  be 
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selected  to  yield  correlations  with  empirically  established  "knockdown" 
factor  versus  effective-radius-to- thickness  curves  as  suggested  by,  for 
example,  Madsen  and  Hoff29  and  Almroth  et  al.2  However,  in  Reference  2 
arguments  are  presented  for  limiting  the  usefulness  of  the  Koiter  theory 
in  practical  analysis  at  the  present  time. 

Effectively  Thick  Cylinders 

In  the  light  of  the  discussion  to  this  point,  it  is  obvious  that  the 
presence  of  initial  imperfections  is  the  main  cause  of  the  discrepancy 
between  prediction  and  experiment  for  thin  shells  under  compressive 
loading.  However,  the  effect  of  imperfections  tends  to  diminish  as 
the  radius-to- thickness  ratio  is  decreased.  Stiffened  and  sandwich 
shells  of  practical  dimensions  are  by  nature  in  the  "effectively 
thick"  category.  Whereas  unstiffened  shells  can  have  "knockdown" 
factors  applied  to  the  classical  buckling  loads  approaching  0.20, 
it  has  been  suggested  that  shells  with  closely  spaced,  moderately 
heavy  stiffening  may  be  designed  with  a  "knockdown"  factor  of  0.75 
applied  to  the  classical  buckling  load  as  long  as  the  material  re¬ 
mains  linear  elastic  and  no  local  skin  buckling  precedes  general  in¬ 
stability  (see  Reference  9)-  Further,  the  extensive  experimental  work 
of  Singer  et  al.^*^*^2  has  given  rise  to  the  contention  that  practical 
stiffened  shells  tend  to  behave  reasons’ ly  closely  indeed  to  the  clas¬ 
sical  theory  and  that  a  sufficiently  conservative  reduction  factor  for 
linear  elastic  material  behavior  lies  at  about  the  0.75  level.  Also, 
Hoff  in  Reference  U  has  commented  to  the  effect  that  the  classical 
theory  for  orthotropic  shells  is  probably  sufficiently  accurate  for 
engineering  calculations  when  the  ef fective-radius-to-thickness  ratio 
of  the  stiffened  cylinder  is  equal  to  or  less  than  200  . 

MAXIMUM  STRENGTH  CONSIDERATIONS 


To  the  designer  concerned  with  weight  problems  in  shell  structures, 
there  is  little  comfort  in  a  status  quo  which  leaves  the  satisfactory 
predictions  of  buckling  loads  under  compression  unresolved  theoretically 
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and  no  other  choices  than  either  the  use  of  semiempirical  criteria 
which  are  admittedly  still  conservative  or  design  by  test,  a  costly 
if  not  prohibitive  procedure.  Since  the  cylindrical  shell  buckling 
problem  has  been  treated  ir  at  least  1600  references  and  still  poses 
a  challenge  even  in  the  homogeneous,  isotropic  case,  it  seems  quite 
reasonable  to  consider  another  facet  of  the  physical  behavior  at 
buckling  which  has  been  almost  deliberately  avoided  for  as  long  as 
experiments  have  failed  to  demonstrate  the  predictions  of  the  clas¬ 
sical  theory.  This  facet  has  to  do  with  the  role  that  nonlinear 
material  behavior  may  play  in  the  shell  buckling  process.  The 
buckling  load  of  a  compressed  shell  is  in  reality  the  maximum  or 
failure  load.  To  deal  with  the  failure  mechanism  on  the  basis  of 
linear  elastic  material  behavior,  especially  in  "effectively  thick" 
shells  of  the  stiffened  or  sandwich  type  which  if  devoid  of  initial 
imperfections  would  undergo  buckling  without  yielding,  is  questionable. 
Perhaps  preoccupation  with  postbuckling  behavior  of  very  thin  shells 
at  large  deflections  and  correlation  of  minimum  postbuckling  equilibrium 
stresses  with  the  low  buckling  stresses  of  actual  cylinders  turned  atten¬ 
tion  away  from  the  maximum  strength  problem.  Or  perhaps  preoccupation 
with  imperfection  sensitivity  of  compressed  shells  kept  attention  cen¬ 
tered  qualitatively  in  this  area.  Imperfect  shells  buckling  in  the 
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yield  region  have  been  discussed  by  Batterman  and  Hutchinson;  how¬ 
ever,  this  type  of  problem  concerns  thick  shells  and  cannot  be  considered 
representative  of  the  aircraft  shell  design  problem.  Remarkably,  Donnell^ 
in  his  earliest  work  on  the  buckling  of  imperfect,  thin  shells  in  com¬ 
pression  developed  a  design  formula  which  based  failure  on  the  arbitrary 
criterion  of  the  maximum  stress  at  buckle  nodes  reaching  the  yield  stress 
of  the  material.  Had  this  concept  of  thin  shell  failure  been  pursued 
more  vigorously,  it  is  possible  that  the  succeeding  more  accurate  large- 
deflection  analyses  of  shells  would  not  have  been  restricted  to  linear 
material  behavior  for  so  many  years. 

In  the  course  of  investigating  the  accuracy  of  von  Karman-Donnell-^ 
shell  theory  for  studying  the  axial  compression  postbuckling  problem  of 
thin  shells,  Mayers  and  Rehfield^  and  Hoff,  Madsen,  and  Mayers^  in 


7 


complementary  efforts  brought  an  end  to  the  search  for  the  minimum  post- 
buckling  stress  a  thin  shell  can  support.  The  former  concluded  as  well 
that  the  shell  postbuckllng  problem  cou^d  not  be  resolved  satisfactorily 
on  the  basis  of  linear  elastic  material  behavior.  Subsequent  investiga¬ 
tions  by  Mayers  and  Wesenberg^'  ^  of  the  nonlinear  behavior  of  initially 
imperfect,  axially  compressed  unstiffened,  stiffened,  and  sandwich  shells 
in  the  initial  buckling  region  (maximum  strength  problem)  have  yielded 
indisputable  evidence  that  the  maximum  load  can  be  strongly  dependent 
upon  the  combined  effects  of  initial  imperfections  and  the  loss  of  the 
classical  linear  relationship  between  stress  and  strain  given  by  Hooke's 
law.  The  phenomenon  is  pronounced  in  either  closely  stiffened  shells 
with  moderately  heavy  stiffening  or  sandwich  cylinders  ("effectively 
thick"  practical  construction)  using  materials  that  are  described  by 
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the  lower  exponents  in  the  Ramberg-Osgood'5  representation  of  uniaxial 
stress-strain  curves  for  aircraft  structural  metals.  Materials  in  which 
the  uniaxial  stress-strain  curve  exhibits  little  range  of  persisting 
initial  modulus  include  stainless  steel,  nickel,  copper,  brass,  and 
nonferrous  superalloys;  annealed  metals  in  general;  metals  at  elevated 
temperature;  and,  to  a  large  extent,  filamentary  composites  with  high- 
performance  fibers. 

Optimization  Considerations 

Normally,  in  the  design  process  using,  for  example,  the  criteria  of 
Reference  9>  a  stiffened  or  sandwich  cylinder  is  assumed  to  be  geo¬ 
metrically  perfect  and  in  a  state  of  uniform  compressive  stress  at 
the  instant  of  buckling  (bifurcation).  The  possibility  of  buckling 

at  stresses  exceeding  the  proportional  limit  stress  is  accounted  for 
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by  applications  of  a  reduced  modulus  after  Gerard.  Then,  the  pre¬ 
sence  of  initial  imperfections  in  the  shell  is  recognized  by  applying 

a  suitable  "knockdown"  factor  obtained  from  an  empirical  design  curve 
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such  as  that  of  Weingarten  et  al.  It  is  important  to  note,  as  pointed 
out  in  References  10  and  11,  that  nonlinear  material  behavior  is  not  an 
unusual  occurrence  in  imperfect  shells  undergoing  bending  and  initial 
buckling  (maximum  strength)  even  though  the  average  compressive  stress 


remains  in  the  linear  elastic  range.  Also,  it  should  be  borne  in  mind 
that  the  use  of  a  reduced  modulus  in  the  "handbook"  approach  does  not 
necessarily  imply  that  yielding  in  the  design  sense  (0.2$  offset  yield 
stress)  has  taken  place.  Thus,  it  is  still  moot  whether  or  not  the 
commonly  used  design  procedure  based  upon  a  gross  idealization  of  com¬ 
pressed  shell  buckling  behavior  is  realistic.  Presumably,  the  pro¬ 
cedure  is  well  on  the  conservative  side  and,  therefore,  justifiable. 

Now,  where  structural  efficiency  is  concerned,  the  degree  of  conserva¬ 
tism  is  important.  The  justification  for  the  relatively  new  design 
procedure  presented  in  Reference  2  is  that,  although  based  upon  semi- 
empirical  considerations,  in  the  absence  of  satisfactory  theory,  it 
yields,  nevertheless,  the  same  or  less  conservative  results  than  other 
methods  in  use.  For  a  maximum  strength  (failure)  analysis  of  an  ini¬ 
tially  imperfect,  stiffened  shell,  it  would  appear  that  bending,  buckling, 
and  postbuckling  in  a  nonlinear  elastic  sense  prior  to  yielding  is  more 
representative  of  the  actual  physical  behavior  than  either  (l)  the  clas¬ 
sical  idealization  of  a  perfect  shell  with  correction  factors  to  account 
for  both  initial  imperfections  and  nonlinear  material  behavior  applied 
in  a  gross  sense  or  (2)  prediction  of  the  reduced  buckling  load  based 
upon  a  linear  material,  initial  postbuckling  analysis  which  takes  into 
account  the  first-order  effects  of  an  axisymmetric  imperfection  pattern, 
the  amplitude  of  which  is  correlated  with  empirical  curves  of  "knockdown" 
factor  versus  radius-to-ef fective-thickness  ratio  for  axially  compressed 
cylinders.  Which  of  the  approaches  leads  to  a  more  realistic  estimate 
of  compressed  shell  behavior  for  optimization  considerations  and  the 
evaluation  of  test  data  is  the  basis  of  the  study  results  to  follow. 

The  study  itself  is  based  upon  the  theoretical  developments  presented 
in  References  10,  11,  and  40  which,  in  turn,  are  outgrowths  of  the 
investigations  reported  in  References  4l,  42,  and  4j.  In  the  latter 
studies,  compressed  plates  were  of  interest;  Reference  43  presents  the 
essentially  exact  solution  to  the  maximum  strength  of  postbuckled  rec¬ 
tangular  plates  and  places  the  well-known  von  Kannan  effective  width 
44 

formula  in  proper  perspective. 
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As  a  basis  for  a  necessarily  limited  optimization  study,  the  longer  of 
the  two  2024-T3  aluminum  alloy,  integral,  externally,  long* tudinally 
stiffened  cylinder  cross  sections  used  in  the  experimental  investiga- 
tion  of  Card  and  Jones  is  adopted  for  the  reference  case;  it  appears 
in  Figure  1  as  Case  0  .  The  uniaxial  stress-strain  curve  for  2024-T3 
aluminum  alloy  is  shown  in  Figure  2.  Although  no  attempt  appears  to 
have  been  made  in  the  program  of  Reference  45  to  consider  the  problem 
of  cylinder  optimization,  the  stringers  were  apparently  sized  and 
located  to  preclude  local  buckling  in  view  of  the  fact  that  none 
occurred.  The  actual  cylinder  behaved  •'s  if  it  were  “effectively 
thick",  since  the  experimental  result  was  9 2$  of  the  calculated  clas¬ 
sical  value.  This  is  reasonable  as  the  effective  radius- to- thickness 
ratio  for  Case  0  is  only  167  •  As  shown  in  Reference  10,  aluminum 
cylinders  under  compression  begin  to  exhibit  the  effects  of  stress- 
strain  curve  nonlinearity  below  about  ~  200  .  The  test 

results  of  Reference  45  in  the  form  of  cylinder  stress-strain  curves 
bear  out  the  fact  that  the  reference  cylinder  was  in  the  nonlinear 
range  at  maximum  load.  Thus,  with  the  strain  data  indicating  no  local 
buckling  prior  to  general  buckling,  the  discrepancy  between  classical 
theory  and  experiment  cannot  be  attributed  to  initial  imperfections 
and  boundary  conditions  alone. 

Now,  under  the  hypothetical  design  consideration  that  no  yielding 
(0.2$  offset  definition)  shall  occur  at  limit  load,  the  buckling 
stress  in  the  cylinder  can  be  raised  appreciably  even  though  the 
nonlinear  portion  of  the  stress-strain  curve  below  the  yield  stress 
is  involved.  On  the  basis  of  simply  supported  ends,  the  Baruch/S inger- 
NASA^’^  uniform  compression  buckling  criterion  yields  18.6  ksi  for 
Case  0  as  indicated  in  Table  I;  for  Cases  1,  2,  and  3,  the  results  are  seen 
to  be  31.5,  31.4,  and  40. 1  ksi,  respectively.  Each  of  the  latter  three 
cases  is  in  the  nonlinear  elastic  range  (nonlinear  elastic  and  inelastic 
are  used  interchangeably  for  convenience  below  the  yield  stress)  and  a 
common  design  procedure  (see,  for  example,  Reference  9)  is  to  apply  a 
plasticity  correction  factor  or  reduced  modulus  (t)  =  V( EtEs)/E2’)  to 
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Figure  1.  Cylinder  Geometries  of  Optimization  Study 


Figure  2.  Stress-Strain  Curve  for  2024-T3  Aluminum. 


the  result.  The  corrected  buckling  stresses  are  shown  in  Table  I  along 
with  the  buckling  stress  predictions  of  the  maximum  strength  theory  and 
analysis  outlined  in  Appendix  I.  In  each  case,  an  initial  imperfection 
amplitude  s11q  =  0.0001  has  been  assumed  in  order  to  provide  a  loading 
path  into  the  nonlinear  elastic  range.  For  Case  0,  it  may  be  noted  that 
the  Baruch/S inger-NASA *  perfect  cylinder  criterion  and  the  imperfect 
cylinder  maximum  strength  analysis  yield  essentially  the  same  result 
since  the  material  behaves  linearly  and  the  imperfection  amplitude  ^11q 
is  so  very  small.  For  Cases  1,  2,  and  3,  however,  the  critical  stresses 
obtained  from  the  two  markedly  different  approaches  are  in  good  agree¬ 
ment  only  for  q  ~  1  .  That  is,  in  view  of  the  stress  levels  involved, 

the  Baruch/S inger-NASA  values  must  be  corrected  for  inelastic  effects; 
the  resulting  inelastic  buckling  stresses  are  seen  to  be  increasingly 
lower  than  the  maximum- strength  analysis  counterparts  as  the  linear 
elastic  buckling  stress  levels  increase.  The  discrepancy  grows  still 
further,  Case  0  included,  because  the  design  recommendation  of,  say. 
Reference  9  includes  a  reduction  of  the  Baruch/S  inger-NASA  buckling 
stresses  by  a  factor  of  0.75  to  account  for  the  possible  presence  of 
initial  imperfections.  The  factor  of  0.75  is  recommended  for  closely 
spaced,  moderately  heavy  stiffening;  it  becomes  as  small  as  0.20  for 
unstiffened  cylinders  of  relatively  large  radius-to-thickness  ratio. 
Although  it  is  true  that  the  maximum  strength  analysis  has  been  based 
upon  an  imperfection  amplitude  of  Hl^  =  0.0001  ,  the  imperfection 
amplitude  required  to  reduce  each  of  the  maximum  strength  analysis 
results  by  0.75  (in  consonance  with  the  handbook  recommendation)  is 
unrealistic  as  indicated  by  the  ^lln  values  listed  for  each  case 
in  Table  I.  Since  for  unstiffened  thin  cylinders,  values  of  ^11^  >0.1 
would  imply  unsatisfactory  fabrication  techniques,  the  Table  lvalues, 
which  are  4  to  10  times  greater,  show  that  the  0.75  "knockdown"  factor 
is  too  conservative  for  application  to  all  cylinders  with  closely  spaced, 
moderately  heavy  stiffening,  at  least  for  2024-T3  aluminum. 

The  weight  of  each  cylinder  referred  to  Case  0  is  also  shown  in  Table  I. 
It  can  be  seen  that  Case  1  is  the  lightest  (-  2 6ja)  and  carries  almost 
more  stress  than  does  the  reference  cylinder  according  to  the 
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Ba- uch/Singer-NASA  perfect  cylinder,  linear  theory  criterion  corrected 
for  inelasticity.  By  comparison,  the  maximum  strength  analysis  shows 
that  for  =  0.0001  (an  essentially  perfect  cylinder),  Case  1  car¬ 

ries  about  Gfjo  more  stress  at  buckling.  On  the  other  hand,  the  next 
lightest  cylinder  (-  16$),  Case  3 >  supports  86$  more  stress  according 
to  the  Baruch/S inger-NASA  criterion  corrected  for  inelasticity,  but  al¬ 
most  110$  more  stress  according  to  the  maximum  strength  theory.  The 
reason  for  this  must  be  attributed  to  the  plasticity  correction  factor 
r)  and  the  linear  buckling  theory  assumption  of  a  uniform  stress  state 
in  the  cylinder  at  the  instant  of  buckling  (bifurcation).  At  high  stress 
levels  and  prior  to  reaching  the  material  yield  stress,  the  plasticity 
correction  factor  q  can  reflect  very  strongly  the  tangent  modulus  of 
the  stress-strain  curve  and  appreciably  reduce  the  linear  elastic  bifurca¬ 
tion  stress.  It  can  be  noted  from  Table  I  that  the  inelastic  buckling 
stresses  for  Cases  1,  2,  and  3  based  upon  the  linear  theory  buckling 
criterion  are  closer  together  than  in  the  purely  linear-theory  predic¬ 
tions;  this  is  the  result  of  the  tangent  modulus  behavior  at  increasingly 
nonlinear  stress  levels. 

Unlike  the  linear  analysis,  the  maximum  strength  analysis  is  based  upon 
the  concept  of  bending  and  buckling  of  cylinders  with  both  large-deflection 
stress  redistributions  and  localized  nonlinear  stress-strain  behavior  oc¬ 
curring  prior  to  and  during  the  buckling  process.  This  behavior  appears 
to  be  more  representative  of  the  failure  mechanism  in  compressed  shells 
and,  at  least  from  the  limited  study  results  reflected  in  Table. I,  it 
shows  that  there  is  a  strong  question  of  conservatism  in  a  common  design 
procedure  of  the  type  recommended  in  Reference  9*  The  conservatism  will 
be  amplified  by  considerations  of  imperfections,  since  the  0.75  "knock¬ 
down"  factor  has  been  difficult  to  justify  by  the  maximum  strength  ap¬ 
proach.  Of  course,  for  designs  operating  in  the  linear  elastic  range, 
for  designs  with  light  stiffening,  and  for  designs  in  which  local  buckling 
is  permitted  to  occur  before  general  buckling,  the  conservative  approach 
may  be  justified.  However,  the  discussion  here  concerns  modern  aircraft 
applications  and,  thus,  the  type  of  structure  implied  is  that  described 
by  closely  spaced,  moderately  heavy  stiffening  which  neither  yields  at 
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limit  load,  buckles  locally  prior  to  general  buckling,  nor  fails  in 
catastrophic  fashion. 

From  Table  I  and  Figure  1,  it  can  be  seen  that  with  a  yield  stress  of 
about  48  ksi  for  2024-T3  aluminum  alloy,  the  combination  of  60  stringers 
and  39  rings  with  slightly  increased  basic  shell-wall  thickness  is  a  great 
deal  more  efficient  than  a  60-stringer  cylinder  reflecting  heavier  stringers 
and  a  lighter  basic  skin.  The  result  is  consistent  with  previous  findings 
based  on  purely  linear  elastic  theory  (see,  for  example,  Reference  16) 
that,  for  a  given  set  of  cylinder  properties,  external  stringers  and  ex¬ 
ternal  rings  are  the  most  efficient  combination  for  compressive  loading 
cases.  Also,  in  Reference  40,  Meller  and  Mayers,  on  the  basis  of  re¬ 
sults  of  a  linear  elastic,  large-deflection  analysis  of  imperfect  cyl¬ 
inders,  show  the  importance  of  wavelength  control  in  the  snapthrough 
buckling  of  unstiffened  compressed  shells  and  the  beneficial  role  that 
circumferential  stiffening  can  play  with  respect  to  critical  wavelengths 
and,  in  turn,  buckling  strength.  As  mentioned  earlier,  the  postbuckled 
state  of  failed  shells  is  important;  whereas  catastrophic  collapse  is 
typical  of  unstiffened  thin  cylinders  under  compression,  the  postbuckled 
state  would  appear  to  be  less  severe  when  closely  spaced  ring  stiffening 
is  used  in  combination  with  stringer  stiffening.  Photographs  of  buckled 
unstiffened,  stringer-stiffened  and  ring-stiffened,  and  ring-  and  stringer- 
stiffened  shells  are  presented  in  References  46,  32,  and  47,  respectively. 

As  a  point  of  interest,  in  Reference  32  it  is  mentioned  that  "the  lightly 
ring-stiffened  and  the  stringer-stiffened  shells  all  buckled  with  a  loud 
snap  whereas  the  heavily  ring-stiffened  shells  exhibited  a  more  gentle 
type  of  buckling".  Unfortunately,  according  to  Peterson,*  there  is  no 
test  data  in  existence  for  buckling  under  compressive  loadings  of 
stringer-  and  ring-stiffened  shells  in  which  local  buckling  does  not 
precede  general  instability.  Further,  according  to  Peterson  again,* 
testing  through  1969  has  not  disclosed  any  type  of  stiffened  cylinder 
construction  that  fails  in  general  instability  under  compression  and  is 
immune  to  low  failing  loads.  Perhaps  the  use  of  "effectively  thick" 
stiffened  cylinders  and  maximum  strength  theory  can  be  useful  in 
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assisting  in  the  resolution  of  this  very  practical  design  optimization 
problem. 

Credibility  of  Maximum  Strength  Analysis 

For  practical  s  ells  of  closely  stiffened  or  sandwich  construction  under 
compressive  loading,  the  effective  radius-to-thickness  ratio  lies  be¬ 
tween  100  and  300.  Below  100,  plastic  buckling  (with  yielding)  is  in¬ 
volved  and  is  of  little  interest  to  the  aircraft  shell  design  problem. 
Above  about  3^0,  the  effects  of  initial  imperfections  and  local  buckling 
lower  the  general  buckling  strength  in  increasingly  deleterious  fashion 
and,  hence,  exact  too  great  a  weight  penalty  from  the  basic  load-carrying 
structure;  this  weight  penalty,  in  turn,  is  amplified  through  the  over¬ 
all  vehicle.  The  critical  situation  being  faced  with  regard  to  basic 
structural  weight  is  evidenced  by  the  continued  reference  for  at  least 
a  decade  to  advanced  composite  materials  as  the  panacea  for  the  struc¬ 
tural  efficiency  problem. 

For  the  "effectively  thick"  cylinders  in  the  effective  radius-to-thickness 
range  100-300,  the  imperfections  tend  to  minimize  their  influence  and  the 
failure  mechanism  involves  both  the  inherent  sensitivity  of  thin  shells 
under  compression  to  initial  imperfections  and  the  nonlinear  behavior  of 
the  material  (see,  for  example,  Hutchinson  and  Koiter  and  Mayers  and 
Wesenberg^’  ,  This  behavior  is  quite  different  from  that  of  flat 
rectangular  plates  under  compression  where  the  failure  mechanism  de¬ 
pends  on  the  phenomenon  of  buckle  wavelength  change  and  the  nonlinear 
action  of  the  material  (see  References  48  and  43);  imperfection  sensi¬ 
tivity,  which  is  important  in  the  bifurcation  region,  has  little  effect 
on  the  maximum  strength  of  postbuckled  plates.  In  "effectively  thick" 
shells,  even  with  small  imperfections  present,  bendirg  and  buckling  in 
a  nonlinear-elastic  sense  occur  prior  to  and  at  maximum  load.  Thus, 
neither  plastic  buckling  of  thick  shells  nor  elastic  buckling  of  very 
thin  shells  is  a  physically  valid  approach  for  "effectively  thick" 
shells,  the  generous  use  of  plasticity  correction  factors  and  initial 
imperfection  "knockdown"  factors  notwithstanding.  The  state  of  the  art 
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Is  assessed  very  succinctly  by  Almroth  et  al.  in  the  statement,  "Al¬ 
though  it  is  felt  that  the  method  recommended  here  represents  a  clear 
advantage  over  present  design  practices,  it  is  still  an  interim  solu¬ 
tion  acceptable  only  because  totally  satisfactory  methods  ate  not 
available". 

For  "effectively  thick"  cylinders,  the  discrepancy  between  experiment 
and  theory  can  be  as  much  as  40$.  This,  of  course,  is  far  less  dis¬ 
crepancy  than  that  found  with  unstiffened  cylinders  where  the  etperi- 
mental  results  may  be  8056  lower  than  the  theoretical  predictions.  With 
a  reduction  factor  of  0.75  recommended  in,  say,  Reference  9  for  shells 
with  closely  spaced,  moderately  heavy  stiffening,  and  for  the  effective 
radius-to-thickness  ratio  range  100-300,  it  is  reasonable  to  assume 
that  both  imperfections  and  material  nonlinearity  affect  the  failure 
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load  in  a  maximum  strength  sense  as  opposed  to  the  concept^  of  perfect- 
shell  bifurcation  buckling  with  applications  of  semiempirical  correction 
factors  for  imperfection  sensitivity  and  material  nonlinearity.  The 
relatively  new  design  procedure  of  Reference  2,  which  includes  the 
Koiter  theory  as  one  of  its  bases,  is  restricted  to  linear  material 
behavior. 

Although  Mayers  and  Wesenberg^’^  have  demonstrated  the  combined  effects 
of  imperfections  and  material  nonlinearity  on  shell  maximum  strength, 
no  other  attempt  has  been  made  to  offer  a  possible  reason  other  than 
imperfections  for  the  large  scatter  appearing  in  plots  of  myriad  non- 
dimensionalized  experimental  buckling  stress  results  versus  radius-to- 
thickness  ratio  (see,  for  example,  References  25  and  4M.  Donnell,  ^ 
in  an  amazing  effort  as  far  back  as  193^,  attempted,  among  other  things, 
to  correlate  the  buckling  strengths  of  unstiffened,  imperfect  steel, 
brass,  and  Dural  cylinders  with  a  semiempirical  criterion  based  upon 
material  yielding  of  buckled  sheet;  he  achieved  good  success  in  the 
radius-to-thickness  range  150-1500.  Since  the  Donnell  theoretical 
developments  were  not  rigorous,  it  is  a  source  of  bewilderment  as  to 
the  lack  of  interest  shown  in  Donnell's  intuitive  reasoning  and  results 
in  the  years  following  when  so  much  effort  was  expended  in  attacking 
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rne  thin  shell,  axial  compression  buckling  problem  and  attempting  to 
explain  both  qualitatively  and  quantitatively  the  discrepancy  between 
experiment  and  theory. 

25 

In  1965,  Weingarten  et  al.  proposed  a  design  formula  for  the  buckling 
of  axially  compressed  shells  based  upon  assessment  of  their  own  and  the 
experimental  results  of  many  other  investigators.  The  formula  recognized 
the  dependence  of  buckling  load  reduction  on  the  radius -to -thickness  ratio 
and,  in  essence,  formed  the  lower  bound  to  data  exhibiting  a  very  large 
scatter  band.  The  Weingarten  et  al.  formula  appears  in  Appendix  II.  On 
the  hypothesis  that  some  of  the  scatter  in  *-est  results  can  be  attributed 
to  nonlinear  material  behavior  during  the  buckling  process,  a  major  mod¬ 
ification  of  the  Weingarten  et  al.  formula  is  made  in  Appendix  II  to 
cover  the  scatter  band  obtained  from  the  experimental  data  presented  in 
References  25  and  46  in  the  form  of  nond imensional  buckling  stress  pa¬ 
rameter  versus  radius- to -thickness  ratio.  The  scatter  band  limits  for 
the  combined  data  are  shown  in  Figure  3.  The  lower-bound  curve,  which 

is  slightly  lower  than  that  of  Weingarten  et  al. .  corresponds  approxi¬ 

mately  to  materials  with  a  proportional  limit  strain  of  0.001,  whereas 
the  upper-bound  curve  corresponds  to  materials  with  a  proportional-limit 
strain  of  about  0.007.  Interestingly,  aircraft  metals  reflect  this  range 
of  values  for  proportional-limit  strain  at  room  temperature.  Some  of 
the  steels  are  on  the  lower  bound;  a  material  such  as  titanium  lies  on 
the  upper  bound.  The  trend  of  the  curves  bears  out  the  trend  obtained 

by  Donnell*^  for  his  own  tests  with  brass  and  steel  cylinders  and  for 
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those  of  Lundquist  with  Dural.  For  a  perfect  cylinder,  the  classical 
buckling  strain  is  inversely  proportional  to  the  radius -to- thickness 
ratio;  thus,  as  shown  in  Figure  3,  in  the  lower  radius- to-thickness 
ratio  range  especially,  the  scatter  band  can  be  traversed  for  a  con¬ 
stant  value  of  radius-to-thickness  ratio  by  designating  materials  in 
increasing  values  of  their  proportional-limit  strains.  This  implies 
that  a  perfect  cylinder  with  a  material  proportional-limit  strain  of 
0.0C1  remains  linear  elastic  only  at  a  radius-to-thickness  ratio  above 
600,  whereas  a  material  with  a  proportional- limit  rstrain  of  0.006  remains 
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Figure  3-  Correlation  of  Material-Dependent  Semiempirical  Relation  of  Appendix  II  With  Scatter 
Band  for  Experimental  Axial  Compression  Buckling  Stress  Coefficients. 


linear  elastic  above  a  radius -to -thickness  ratio  of  100.  Naturally, 
actual  cylinlers  are  not  perfect;  therefore,  bending  and  buckling  in 
a  nonlinear-elastic  sense  take  place  with  imperfections  and  nonlinear 
material  behavior  both  contributing  to  the  failure  mechanism  in  im¬ 
perfect  shells.  The  difference  in  material  stress-strain  curves  must 
manifest  itself  and  contribute  to  the  scatter  in  test  results.  This 
conclusion  is  consistent  with  the  maximum  strength  results  obtained 
by  Mayers  and  Wesenberg  in  References  10  and  11  for  materials  such  as 
electroformed  copper  and  nickel,  aluminum,  and  stainless  steel. 

Another  test  of  the  credibility  of  the  maximum  strength  approach  is 

available  through  the  experimental  work  of  Weller  et  al.  on  longi- 

12 

tudinally  stiffened  steel  cylinders.  In  a  series  of  88  tests  con¬ 
ducted  on  5-  and  7-in. -radius  shells  with  essentially  the  same  effective 
thicknesses  in  the  stiffening  direction,  the  investigators  obtained  the 
rather  anomalous  result  that  on  the  average  the  thinner  (l4-in.)  cylin¬ 
ders  were  stronger  in  buckling  than  the  thicker  (10-in.)  cylinders  when 
the  experimental  buckling  stresses  were  ratioed  to  the  buckling  stresses 
predicted  by  the  classical,  eccentrically  stiffened,  orthotropic  shell 
theory. The  buckling-stress  ratios  as  well  as  other  pertinent  de¬ 
tails  are  given  in  Table  II  for  4  of  the  12  cylinders  in  Reference  12 
having  length- to-radius  ratios  equal  to  or  greater  than  unity;  the  cyl¬ 
inder  designations  are  those  of  Reference  12.  Also  included  in  Table  II 
are  the  buckling-stress  ratios  based  upon  the  theoretical  predictions  of 
the  present  maximum  strength  theory  for  the  4  cylinders  using  the  actual 
stress-strain  curves  for  the  materials  involved  shown  in  Figure  4.  Al¬ 
though  the  steel  materials  used  in  the  tests  exhibit  the  identical 
Young's  modulus,  the  nonlinear  portions  of  the  curves  are  significantly 
different.  With  the  weaker  material  assigned  to  the  5-in.  radius  cyl¬ 
inders,  the  results  as  presented  in  Reference  12  tend  to  contradict  all 
previous  experience.  However,  in  the  light  of  the  corrected  buckling- 
stress  ratios  given  in  Table  II  for  the  subject  cylinders,  it  can  be 
seen  that  the  two  7- in.- radius  cylinders  retain  the  same  stress  ratios 
(the  behavior  is  truly  linear  elastic)  whereas  marked  changes  occur  in 
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TABLE  II.  COMPARISON  OF  EXPERIMENTAL  BUCKLING  STRESS  COEFFICIENTS 


the  stress  ratios  for  the  5-fn.- radius  cylinders  due  to  very  definite 
inelastic  behavior.  Indeed,  10-in.  cylinders  SZ-3  and  J-L  reverse  the 
anomalous  trend  appearing  in  Reference  12  with  respect  to  l4-in.  cyl¬ 
inder  24-L;  however,  although  the  gap  is  closed  significantly  with  re¬ 
spect  to  cylinder  13-L,  the  very  heavy  stiffening  of  the  latter  relative 
to  that  of  SZ-3  and  7~L  (see  Table  II)  prevents  an  actual  reversal  of 
the  trend.  Had  cue  assignment  of  material  to  the  10-  and  l4-in.  cyl¬ 
inders  been  reversed  prior  to  fabrication,  the  l4-in.  cylinders  would 
have  appeared  to  be  much  weaker  on  the  average  than  the  10- in.  cylin¬ 
ders.  The  tendency  then,  no  doubt,  would  have  been  to  attribute  the 
discrepancy  entirely  to  imperfection  sensitivity  of  the  "effectively 
thinner"  14-in.  cylinders.  Another  point  of  interest  is  that  the 
question  of  yielding  of  the  material  in  the  5-in.  -radius  cylinders 
is  not  involved.  The  average  compressive  stresses  are  well  below 
the  O.lf,  offset  yield  stress  of  the  material  for  the  5-in. -radius 
cylinders.  The  maximum  local  stresses  at  buckling  are  somewhat 
higher,  although  still  below  the  yield  stress.  Postbuckling  local 
maximum  stresses  are  another  problem  when  catastrophic  snapthrough 
occurs;  such  considerations  are  beyond  the  scope  of  the  present  dis¬ 
cussion.  However,  they  are  discussed  in  Reference  $6. 

The  discussion  and  results  presented  in  this  section  are  believed  to 
lend  credence  to  the  concept  of  a  maximum  strength  analysis  of  com¬ 
pressed  cylinders.  The  results  are  consistent  with  the  conclusions 
of  Mayers  and  Wesenberg^  *  that  the  effects  of  material  nonlinearity 
together  with  those  of  initial  imperfections  control  the  buckling 
strength  of  compressed  cylinders.  Material  nonlinearity  plays  the 
more  significant  role  quantitatively  for  "effectively  thick"  cylinders 
(closely  stiffened  and  sandwich  construction),  especially  when  the 
material  exhibits  a  large  ratio  of  yield  stress  to  proportional  limit 
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stress  (or  a  low  value  of  the  exponent  appearing  in  the  Ramberg-Osgood'3 
representation  for  uniaxial  stress-strain  curves  of  conventional  air¬ 
craft  metals). 
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The  experimental  results  of  Weller  et  al.  in  Reference  12  when  inter¬ 
preted  in  the  light  of  the  maximum  strength  concept  of  compressed  cyl¬ 
inder  failure  advanced  in  References  10  and  11  as  well  as  in  the 
present  effort  demonstrate  the  sensitivity  of  the  buckling  process 
to  both  material  nonlinearity  and  initial  imperfections.  The  former 
phenomenon  is  vividly  illustrated  for  the  two  5“in. -radius  cylinders 
taken  from  Reference  12  by  noting  the  experimental  buckling  stress 
levels  given  in  Table  II  with  respect  to  the  stress-strain  curves  of 
Figure  4.  Offhand,  it  would  appear  that  the  buckling  stresses  (or 
average  axial  stresses  at  buckling)  carried  by  the  10-in.  cylinders 
are  essentially  in  the  linear  elastic  range.  However,  the  maximum 
strength  analysis  for  imperfect  shells  involves  both  bending  and 
buckling  (rather  than  the  classical  concept  of  bifurcation  from  a 
perfect  shell  membrane  state  of  uniform  stress)  so  that  local  effec¬ 
tive  stress  levels  (membrane  plus  bending)  can  be  well  into  the  non¬ 
linear  range  of  the  stress-strain  curve.  Failure  is  an  integrated 
effect  for  the  shell  buckling  problem;  thus,  buckling  based  on  the 
maximum  strength  analysis  for  the  10-in.  cylinders,  although  outwardly 
appearing  to  be  essentially  linear  elastic,  actually  is  significantly 
inelastic.  It  is  obvious  that  the  l4-in.  cylinders  remain  linear 
elastic  during  the  bending  and  buckling  process  because  of  the  lower 
stress  levels  involved  and  the  increased  linear  portion  of  the  ma¬ 
terial  stress-strain  curve  for  these  cylinders.  One  more  point  of 
note  is  that  the  nonlinear  material  behavior  begins  to  make  its  in¬ 
fluence  felt  in  the  vicinity  of  the  proportional  limit.  The  fact 
that  average  stress  levels  may  be  well  below  the  0.2$  offset  yield 
stress  for  a  material  does  not  justify  treating  the  buckling  problem 
on  the  basis  of  linear  elastic  material  behavior  up  to  the  yield  point 
stress  level.  The  assumption  of  a  perfectly  elasto-plastic  material 
(with  ~r  without  strain  hardening)  is  a  case  in  point.  Obviously  there 
is  moLe  co  the  perplexing  behavior  of  axially  compressed  thin  shells 
than,  for  example,  the  effects  of  imperfection  sensitivity,  prebuckling 
deformations,  and  boundary  conditions.  It  is  felt  that  more  attention 
will  have  to  be  given  to  material  nonlinearity  effects,  especially  if 


the  optimization  with  respect  to  weight  and  cost  of  stiffness-critical 
shell  structures  is  to  become  a  reality  for  conventional,  sandwich, 
and  advanced  composite  materials  at  both  room  and  elevated  temperatures. 
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A  maximum  strength  analysis  of  initially  imperfect,  axially  compressed, 
eccentrically  stiffened  cylinders  based  upon  a  modification  of  Reissner's 
variational  principle,  von  Karman- Donne 11  shell  theory,  and  a  deformation 
theory  of  plasticity  is  undertaken  to  optimize  in  gross  fashion  a  2024-T3 
aluminum  alloy  configuration  subject  to  the  condition  that  no  yielding 
take  place  in  the  design  sense.  The  results  are  compared  with  the  pre¬ 
dictions  of  linear  theory  for  perfect  shells  with  applied  buckling  load 
reduction  factors  to  account  for  imperfections  and  inelastic  behavior. 

For  the  optimized  cylinder,  combining  stringer-  and  ring-stiffening,  the 
maximum  strength  analysis  is  shown  to  be  less  conservative. 

To  lend  credence  to  the  maximum  strength  approach,  the  technique  is  used 
to  correlate  the  results  obtained  for  several  10-  and  l4-in.  diameter 
cylinders  constructed  of  steel  with  different  proportional  limits,  re¬ 
spectively.  The  maximum  strength  results  explain  the  apparently  anomalous 
behavior  displayed  by  the  two  sets  of  cylinders  when  the  experimental 
buckling  stresses  are  ratioed  to  the  elastic  theory  predictions  of  clas¬ 
sical  theory  for  eccentric  stiffening.  It  is  demonstrated  that  material 
nonlinearity  is  an  important  factor  in  buckling  load  calculations  of  im¬ 
perfect  shells  at  average  compressive  stresses  well  below  the  material 
yield  stress,  especially  in  materials  exhibiting  lower  exponents  in  the 
Ramberg-Osgood  representation  of  the  uniaxial  stress-strain  curve  for 
aircraft  metals.  This  factor  is  of  consequence  for  materials  at  elevated 
temperatures  and  for  composites  using  high-strength  filaments  such  as 
boron. 
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The  effect  of  material  nonlinearity  in  the  effective  radius-to-thickness 
range  characterizing  closely  spaced,  moderately  heavy  stiffened  and  sand¬ 
wich  shells  suggests  that  it  may  serve  to  explain  some  of  the  large  scatter 
appearing  in  experimental  data  which  is  presently  attributed  mainly  to  in¬ 
itial  imperfections.  A  simple,  material-dependent  relation  for  buckling 
coefficients  is  offered  to  show  that  the  scatter  band  encompassing  myriad 
teat  results  can  be  spanned  at  any  radius-to-thickness  ratio  below  1500 
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as  a  function  of  proportional-limit  strains  in  the  range  0.001  through 
0.007.  Interestingly,  this  range  covers  the  proportional  limit  strains 
for  conventional  metals.  The  qualitative  sensitivity  of  the  buckling 
coefficient  relationship  to  material  nonlinearity  cannot  be  overlooked 
when  attempting  to  explain  the  discrepancy  between  test  and  theory  for 
axially  compressed  shells  in  the  "effectively  thick"  range  (practical 
stiffened  and  sandwich  shells).  The  work  of  Donnell  in  1934  in  this 
respect,  although  open  to  question  with  regard  to  rigor  of  analysis 
and  adopted  failure  criterion,  was  apparently  too  soon  overlooked  in 
the  massive  theoretical  attack  on  the  small  and  large  deflection  stability 
problems  of  perfectly  elastic  cylinders  which  has  taken  place  in  the  past 
30  years. 

In  view  of  the  desire  for  increased  structural  efficiency,  the  concept 
of  bending  and  buckling  in  the  inelastic  range  well  below  the  0.2$  off¬ 
set  yield  stress  should  be  given  some  consideration  for  failure  load 
determination  when  "effectively  thick"  shells  (closely  spaced,  moderately 
heavy  stiffening  and  sandwich  cylinders)  are  involved.  Optimized  stringer- 
and  ring-stiffened  shells  should  not  exhibit  catastrophic  snapthrough  col¬ 
lapse  at  limit  load  in  the  sense  associated  with  very  thin,  unstiffened 
cylinders.  Further,  local  buckling  of  the  skin  between  stiffening  ele¬ 
ments  should  not  occur  prior  to  general  buckling.  Unfortunately,  there 
is  no  test  data  for  stringer-  and  ring-stiffened  shells  in  which  there 
is  an  absence  of  a  large  discrepancy  between  theory  and  experiment;  thus, 
this  gap  in  experimental  data  should  be  closed.  Apparently,  the  effect 
of  local  skin  buckling  prior  to  general  buckling  reduces  stiffness  for 
the  latter  more  than  would  be  desired  in  efficient  design.  Additionally, 
there  is  strong  evidence  in  the  present  work  to  suggest  that  the  axial 
compression  test  data  scatter,  presently  attributed  mainly  to  initial 
imperfections,  may  be  appreciably  sensitive  to  material  nonlinearity  at 
buckling  in  the  radius-to-thickness  range  of  practical  interest.  A 
modest  experimental  program  is  recommended  for  constant  radius-to- 
thickness  ratio  cylinders  of  several  materials  exhibiting  proportional 
limit  strains  from  0.001  (steel)  to  0.007  (titanium)  wherein  the 
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fabrication  process  is  identical.  In  this  manner,  some  assessment  of 
the  magnitude  of  scatter  due  to  material  nonlinearity  effects  can  be 
obtained,  other  than  that  implied  from  the  stringer-stif fened-cylinder 
test  program  of  Reference  12  and  the  presentation  contained  in  Appen¬ 
dix  II, 


f 
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APPENDIX  I 


FOUNDATION  OF  MAXIMUM  STRENGTH  ANALYSIS  FOR  INITIALLY 
IMPERFECT,  ECCENTRICALLY  STIFFENED,  CIRCULAR 
CYLINDRICAL  SHELLS  UNDER  AXIAL  COMPRESSION 


For  completeness,  the  cyl  ndrical  shell  maximum  strength  analysis  de¬ 
veloped  in  Reference  10  and  extended  in  Reference  11  is  presented  here 
in  concise  form.  It  is  based  upon  the  use  of  a  modified  version  of  the 
Relssner  variational  principle,  the  von  Karman-Donnell  strain-displacement 
relations,  and  a  deformation  theory  of  plasticity.  Because  the  deflec¬ 
tions  and  rotations  at  maximum  load  are  sufficiently  small,  as  determined 
by  Mayers  and  Wesenberg  in  Reference  10,  the  use  of  the  von  Karman-Donnell 
theory  is  justified.  The  effects  of  unloading  can  be  neglected; ** 
hence,  the  analysis  reduces  to  one  of  nonlinear  elasticity.  To  avoid 
the  complexity  of  integrating  nonlinear  stress  distributions  through 
the  thickness,  a  two-element  cross  section  with  an  intermediate  core 
rigid  in  shear  is  used  to  replace  that  of  the  stiffened  shell  (see 
Figure  5).  The  stiffening  elements  are  sufficiently  closely  spaced 
so  that  they  may  be  considered  Bmeared  over  the  shell  surface,  thus 

giving  rise  in  the  model  to  implied  facing  separations  h  ,  h  ,  and 

x  y 

h  .  The  eccentricty  effects  of  any  stiffening  elements  are  accounted 

for  approximately  in  further  definitions  of  the  quantities  h  ,  h  , 

x  y 

and  h 

xy 


The  Relssner  functional  for  the  orthotropic  cylinder  model  under  con¬ 
trolled  end  shortening,  modified  to  include  incompressible  plasticity 
of  the  material  (all  Poisson  ratios  are  taken  equal  to  £),  is 
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The  orthotropic  shell  constants  take  the  following  form  for  eccentric 
ring-  and  stringer-stiffening. 
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Again,  it  should  be  noted  that  in  the  actual  computations  of  maximum 
strength,  all  Poisson  ratios  are  taken  equal  to  \  . 
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The  parameters  K  and  N  appear  in  the  Ramberg-Osgood^  representa¬ 
tions  of  the  uniaxial  stress-strain  curve  for  aircraft  metals;  namely, 


■(:)•■(;) 


(38) 
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As  shown  in  References  10  and  11,  the  uniaxial  stress-strain  curve 
enters  the  maximum  strength  analysis  through  the  relationship 


The  displacement  distributions  assumed  for  the  imperfect,  uniformly 
shortened  shell  are 
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where  e  ,  |11q  ,  |20o  ,  Xx  ,  and  Xy  are  prescribed,  and  S20q  =  Ull0A) 
for  the  calculations. 

With  no  arbitrary  parameters  appearing  in  the  middle-surface  displace¬ 
ments  u  and  v  ,  middle-surface  equilibrium  is  automatically  satisfied 
by  the  following  set  of  middle-surface  stresses: 
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t 


and 


The  bending  stresses  are  taken  in  the  form  of  the  curvatures  as 


The  modified  Reissner  functional  has  now  been  cast  in  terms  of  IT  free 
parameters;  Sqq  *8  obtainable  from  enforcement  of  the  periodicity 
condition  on  the  v  displacement.  The  vanishing  of  the  first  vari¬ 
ation  of  the  Reissner  functional  leads  to  IT  equations  in  the  IT 
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unknown*.  For  a  given  material  (K,N)  ,  geometry  (R/t  )  >  an<*  Prc" 

•  cribed  value*  of  e  >  Ihq  >  M  =  ^xAy  »  anc*  T  a  n  t^/R  ,  a  load- 
shortening  relationship  in  the  form  crR/Etx  versus  eR/tx  can  be 
obtained;  the  load-shortening  curve  possesses  a  distinct  maximum 
except  for  cylinders  in  which  the  maximum  amplitudes  of  the  initial 
Imperfection  pattern  are  unreasonably  large.  For  the  linear  elastic 
case,  the  system  can  be  reduced  to  four  equations  in  the  four  unknowns 

^11  ’  ^ 20  ’  anc* 
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APPENDIX  II 


STRESS-STRAIN  CURVE  NONLINEARITY:  A  POSSIBLE 
CAUSE  OF  SCATTER  IN  TEST  RESULTS  FOR 
BUCKLINC  OF  CYLINDERS  IN  COMPRESSION 


The  large  amount  of  scatter  that  exists  In  plots  of  experimental 
buckling  atress  coefficient  crR/Et  versus  radius-to-thickness  ratio 
R/t  is  illustrated  in  Figure  3*  The  scatter  band  shown  encompasses 
the  extensive  date,  accumulated  in  References  25  and  46.  The  lower 
bound  to  the  data  of  Reference  25  is  given  by  the  empirical  relation 
developed  by  Weingarten  et  al.  in  Reference  25;  namely, 

C  =  25  a  0.606  -  0.546  j  1  -  exp 


-*VT 


(50) 


where  the  first  term  on  the  right-hand  side  represents  the  classical 
buckling  stress  coefficient  1/V^  -v*" )'  for  Poisson's  ratio 

v  =  0.30  . 

A  modified  relation  is  suggested  for  present  purposes  as 


0.606  -  0.546  (  1  -  exp 


(51) 


Equation  (51)  corresponds  closely  to  the  scatter  band  of  Figure  3  at 

the  limits  e  ,  0.001  and  e  ,  -  0.007  •  Aircraft  structural 

pi  pi 

metals  in  general  possess  proportional  limit  strains  in  the  range 

0.001  <  e  .  <  0.007  •  Equation  (51)  can  be  recast  to  include  the 

pl  13 

parameters  used  by  Donnell  in  1934  in  the  form 


3 


o-R 

C  =  —  =  0.606  -  0.546  jl  -  exp 

Et 
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For  reference,  Donnell '«  dcalgn  equation  taken  from  Reference  13  la 

(53) 

Et  ll  +  0.004  ( E/cry )  I 

It  can  be  observed  that  Equation  (52)  includes  one  more  parameter, 

(T y/°pj  »  than  does  Equation  (53 )•  When  the  ratio  cry/c approaches 
unity,  only  the  parameters  perceived  by  Donnell  appear.  For  c ^/tr  j 
approaching  unity,  however,  the  equations  cannot  be  compared  in  gen¬ 
eral,  since  they  plot  with  opposite  curvatures.  Normally,  for  aircraft 
materials,  the  ratio  ^  is  about  4/3  or  greater. 

For  several  values  of  ,  Equation  (51)  is  plotted  in  Figure  3*  It 

can  be  seen  that  the  scatter  band  is  traversed  at  any  R/t  ratio  as  a 
function  of  material  nonlinearity.  If  it  is  assumed  that  Donnell' s 
brasa  cylinders  possessed  an  e  ^  of  0.0015  (it  i*  known  that  the 
yield  stress  was  only  about  27  kst)  contrasted  with  0.003  for 

the  Dural  cylinders,  then  the  spread  between  the  two  curves  correlates 
very  well  with  the  actual  average  behavior  of  the  21  brass  and  45  Dural 
test  cylinders  for  the  range  300  <  R/t  <  1500  (see  Reference  13). 

The  strong  influence  of  material  nonlinearity  indicated  by  the  develop¬ 
ments  in  this  appendix  cannot  be  overlooked  in  attempts  to  explain  the 
large  scatter  present  in  test  data  for  axial  compression  buckling  of 
thin  shells.  The  scatter  heretofore  has  been  attributed  principally 
to  initial  imperfections  through  theoretical  analyses  involving  linear 
material  behavior.  It  appears,  however,  that  the  question  has  not  been 
answered  as  to  just  how  much  influence  on  scatter  may  result  from  ma¬ 
terial  nonlinearity.  After  all,  as  Donnell  suggested  almost  40  years 
ago,  the  buckling  (or  failing)  stress  of  an  imperfect  thin  cylinder 
in  compression  should  be  a  function  of  the  yield  stress.  The  conten¬ 
tion  herein,  as  well  as  in  the  forerunning  studies  by  Mayers  and 
Wesenberg, **  is  that  buckling  (or  failure)  of  an  imperfect  com¬ 
pressed  cylinder  is  a  function  of  material  nonlinearity  at  stresses 


46 


well  below  the  yield  atreaa,  especially  In  the  range  of  "effectively 
thick"  cylindera  typical  of  cloaely  atiffened  and  aandwich  conatruction. 
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